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Abstract Following Einstein's definition of Lagrangian density and gravita- 
tional field energy density [1-3], Tolman derived a general formula for the total 
matter plus gravitational field energy (Po) of an arbitrary system[4-6]. For 
a static isolated system, in quasi-Cartesian coordinates, this formula leads 
to the well known result P = J ^/ Z7 g(T$ - - T 2 2 - T 3 3 ) d 3 x, where g is 
the determinant of the metric tensor and T h a is the energy momentum tensor 
of the matter. Though in the literature, this is known as "Tolman Mass", it 
must be realized that this is essentially "Einstein Mass" because the underly- 
ing pseudo-tensor here is due to Einstein. In fact, Landau -Lifshitz obtained 
the same expression for the "inertial mass" of a static isolated system with- 
out using any pseudo-tensor at all and which points to physical significance 
and correctness of Einstein Mass [7]! For the first time we apply this general 
formula to find an expression for Po for the Friedmann- Robertson -Walker 
(FRW) metric by using the same quasi-Cartesian basis. As we analyze this 
new result, physically, a spatially flat model having no cosmological constant 
is suggested. Eventually it is seen that conservation of Po is honoured only in 
the a static limit. 
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1 Introduction 



Einstein was the first relativist to point out that in curved spacetime matter 
energy momentum itself is not conserved [1-3]. On the other hand, what is 
conserved is the joint energy momentum of matter and the gravitational field 
produced by matter itself. This is somewhat like the fact that for a system 
of interacting charges (matter), it is the sum of the matter energy momen- 
tum and the electromagnetic field energy momentum produced by the charges 
themselves which is conserved. However, in the presence of gravitation, be- 
cause of curved spacetime what is eventually conserved is a sort of geometry 
weighted energy momentum density which is now called "Energy Momentum 
Complex" (EMC) [6, 8]: 

The Einstein EMC is given by [3-6] 



9 b a = y ^(T b + t b a ) 



(2) 



Here t b a is the gravitation field energy momentum (pseudo) tensor constructed 
from metric tensor g a b and dg a b/dx c . The "energy complex" derived from 
Einstein EMC obviously is 



(3) 



where the subscript referes to time coordinate. The Einstein pseudotensor 
is often referred to as the "canonical" pseudotensor because it is derived by 
using the general formula for the energy momentum of a classical field with a 
Lagrangian density and field variables which could be tensors of any rank. In 
the given case, the field variables are the components of g a t. Obviously, because 
of its the nontensor nature, the local energy density of the field does not have 
a covariant significance. Indeed, this formulation of the consevation law got 
criticized because, in spherical coordinates, empty Minkowski spacetime would 
appear to possess finite energy density. Einstein tried to answer such criticism 
in two ways [9]: 

(i) He showed that though t\ is not a true tensor, the total energy momen- 
tum 

Pa = J el d\ (4) 

nevertheless behaves as a free vector (i.e., a vector not defined at a particular 
spacetime point) under linear transforformations for an asymptotically flat 
spacetime (AFST). 

(ii) He also pointed out that, one must evaluate t b a and 9 b a in quasi-Cartesian 
coordinates. When this prescription is followed, an empty Minkowski space- 
time indeed yields nil value for both these quantities (however see later dis- 
cussion). 
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Since then, all the authors have recognized the fact that for a meaning- 
ful physical answer, the Einstein pseudo tensor must be evaluated in quasi- 
Cartesian coordinates [6, 8]. A case in point would the question of energy mo- 
mentum flux associated with cylindrical gravitational waves (GW), probably 
the only known exact solution of GW. For instance, initially Rosen computed 
the energy momentum flux for this problem using 9 h a (Einstein) in symme- 
try adapted cylindrical polar coordinates and found nil flux [10]. This led him 
to conclude that GWs may not be carrying any energy momentum flux and 
hence, in a sense, pseudo or fictitious quantity. But following the suggestion 
of S. Virbhadra, when Rosen & Virbharda recomputed the energy flux using 
quasi-Cartesian coordinates, they found finite value of flux[\\]. 

Though Tolman too acknowledged the fact that there is an ambiguity in 
the localization of gravitational energy momentum, he [4] claimed that 

"It can be shown, nevertheless, that the equations have the necessary fun- 
damental property of being true in all sets of coordinates, and a completely 
satisfactory justification of the formulation was finally given by Einstein in 
1918[9]." 

In the context of the criticism of the energy momentum conservation Eq.(l) 
which does not involve true tensors, Tolman[4] went on adding that 

"It should be remarked, however, that the appropriate criterion for the 
fundamental significance of equations should not be that they are written in 
tensorial form but that they are written in covariant form so as to be true in 
all sets of coordinates. All tensor equations are indeed covariant equations, but 
this does not exclude the possibility of covariant equations, such as the above 
(i.e., Eq.[l]), which are not tensorial. To assume the contrary would be the 
fallacy of the Dormouse in Alice in Wonderland, who said: "I breathe when I 
sleep" is the same as "I sleep when I breathe." 

However, the fact remains that even if Eq.(l) may look covariant, Eq.(2) 
is certainly not so. It may be also worthwhile to repeat the well known fact 
that, since energy momentum conservation after all appears as a divergence, 
the choice of the pseudo tensor or EMC is by no means unique. In general, the 
EMC may have a form [8] 

6 b a (general) = (- g r^(T b a +t» a ) (5) 

where n is a positive integer. An interested reader may look into the forms of 
some well known pseudo-tensors in the literature[6,8]. This freedom of choos- 
ing various EMCs could be something like the gauge freedom in choosing vector 
potentials in electrodynamics. In particular, the M0ller pseuotensor, by con- 
struction, is a true tensor and thus yields a coordinate independent covariant 
description [12]. However, this is also the only pseudo-tensor which contains 
2nd derivatives of g a b. 

In general, in view of the non-covariant nature of the pseudo-tensor based 
energy-momentum localization, the coordinate independent quasilocal approach 
may be conceptually more important [13]. In this approach, quasilocal energy 
momentum is obtained from an appropriate Hamiltonian. However, by no 
means, this approach solves the ambiguity with respect to the assignment 
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of gravitational energy momentum because there could be infinite choices of 
Hamiltonian boundary conditions. Further, Chang & Nester pointed out that 
the quasilocal approach is intrinsically linked to the pseudotensor approach be- 
cause each pseudotensor seems to be associated with a legitimate Hamiltonian 
boundary condition[14]. Thus the importance of the pseuotensor approach may 
not be undermined on the plea that it does not offer a covariant description. 
As to the criticism, that, all pscudo tensor densities should vanish if evalu- 
ated in Riemannian normal coordinates where first derivatives of g ab vanish, 
one may appreciate that such a vanishing could be a necessary consequence of 
principle of equivalence on which general relativity is based. In particular, as 
pointed out by Xulu[6] and as to be pointed out by us later, Einstein pseudo 
tensor, when evaluated in quasi-Cartesian coordintes, may indeed yield a good 
physical description. 

The field part in the Einstein EMC is given by [3-6] 

where A is the cosmological constant. Note that Tolman did include A in this 
context; for instance see Eq.(87.11) of Ref.[5]. Here the Lagrangian density is 
given by 

l = V~9 g ab (r p ab r« q r* r&) (7) 

in terms of the Christoffel symbols 

r& = \a ad {dcgdb + d b9cd - d d9bc ) (8) 

Note, here, d a denotes differentiation with the coordinate x a . Also, the symbol 



ab 



Recall that the Eqs.(l) and (2) look similar to flat spacetime electromagnetic 
counterpart [5-7]: 

^[T b a (m)+t(f)]=0 (10) 

where T b (m) is due energy momentum tensor of matter (i.e., charge) and t b (f) 
is the same for the electromagnetic field. If we would define 

Pa= J^(T^+t a )d 3 X (11) 

where denotes time coordinate and the integration spans entire 3-space. 
Tolman[4,5] showed that for a finite system with a boundary in AFST 

8P 

^ = Flux of (T aa + t aa ) (12) 
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Thus, following Einstein and Tolman, atleast for isolated objects lying in 
AFST, it seems natural to define P a as the energy- momentum 4- vector of the 
system even though P a is not a true covariant vector[l-6]. Let us now call 



d 3 x (13) 



as the total "so-called" Einstein energy of the system. We use this adjective 
because Pq still depends on coordinates. 

As pointed out by one anonymous referee, in the presence of a A, one 
would not have an AFST, and, in a strict sense, the above notion of "flux" 
and associated energy momentum conservation will break down. However, by 
definition "universe" is the ultimate "isolated" object from which there cannot 
be any outward flux. Then, we may indeed expect energy momentum of the 
universe to be conserved by definition. Further since the FRW metric does 
admit a time like (conformal) Killing vector [15, 16], we may expect properly 
defined total energy to be conserved. 

Few years after Tolman's seminal work in 1931 [4], von Freud[17] showed 
that Einstein EMC could be expressed as the divergence of an appropriate 
superpotential: 

" — 
16tt 



where a comma denotes partial diffentiation and the super-potential is given 

by 

H a b c = -j= g g ap [-g(g bp g cq - g cp g bq )] >q (15) 

One may further note that in terms of this superpotential, Einstein energy 
complex is given by 

e ° = Tir^c (is) 

In this superpotential formalism, however, the physically important quan- 
tity T\ does not appear explicitly. In contrast, as we would see, 6 b a , when 
evaluated directly would contain T* through the Ricci tensor R b a . In the fol- 
lowing we review this formalism and point out that atleast for a static system, 
the result obtained by this generic formalism matches with one which does 
not invoke pseudo potential at all[7\. Finally, we will use this Tolman ansatz 
to find an expression for P for the FRW universe directly by retaining the 
physically important T% term. It may be however be stressed that, the Tol- 
man ansatz does not involve any new psedotensor. On the other hand, it is an 
ansatz to find Pq directly by using Einstein psedotensor. Therefore, though, in 
the literature, the resultant P is often referred to as the "Tolman Mass" , we 
would insist for the term "Einstein Energy" for the same. 
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2 Einstein - Tolman Formula 

Let us recall the Einstein Eq. 

R b a -\Rg b a + Ag b a = -^T b a (17) 

where 

R = 8tt (T ° + Tl + Tl + T|) + 4/1 (18) 

is the Ricci scalar. By combinging the Einstein Eq. and the expression for 
t b a (Eq.[6]), Tolman showed that (see Eq.89.1)[5], one can write the Einstein 
EMC as 

Einstein) = + ±£^=9* - + ±.%L (19) 

where L is related to vacuum Einstein Lagrangian density Lg — yJ—gR, in 
the following way [4-6]: 

L = dc {^^f) ^~ 9R (20) 
By inserting this expression for L, we finally obtain 

6 b a (Einstetn) = -^~V^gR b a - ^g™4ia + ^S b a d c ( g at '^r ) (21) 
ay ' 8tt y 11 16tt a dgl q 16tt a V dgf J v ' 

It should be noted that although A does not explicitly appear in the foregoing 
relation, it remains hidden within the R b term. Thue we obtain 

^Einstein) - - + (,-|L) ,22) 

Since, 

R° = 4w(Tl + Tl + Tl - T °) + A (23) 
we may split 8q into 3 terms: 



where 



9q = A + B + C (24) 



A=\v=g fo-T}-tf-Ti-±} = (-^IlQ), (25) 



16tt 



3siV V V dgf 



(26) 



and 



C = — g ab d ( ] (27) 
16tt V dqf • 
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Since no assumption or precondition has been imposed for the derivation of 
the above formula, it is valid for arbitrary system, whether it has a boundary 
or not, whether it is static or not and whether it is spherical or not. Note the 
physical and geometrical significance of A in view of the presence of and 

in it. Thus for any arbitrary system, spherical or non-spherical, static or 
non-static, having boundary or no-boundary, one may express (see Eq.[23] of 
[4] and Eq.[92.2] of [5]): 

P = j ' 0° d 3 x = J(A + B + C) al 3 x (28) 

where the integration spans over entire 3-space. 



2.1 Finite System with a Boundary 

First application of this formula was made by Tolman[4,5] himself for finding 
the mass-energy of a static finite systen resting in an AFST. Far away from the 
body, the metric is expected to assume the form of an exterior Schwarzschild 
one: 

e v = e~ x = 1 - 2M/R (29) 

where M is an integration constant and R is the circumference coordinate 
(not to be confused with Ricci scalar). Eventually, at very large R, in quasi- 
Catersian coordinates, one should be able to approximate 

als 2 = (1 - 2M/R)dt 2 - (1 + 2M/R)(dx 2 + dy 2 + dz 2 ) (30) 

In this AFST, derivation of Kepler's 3rd law shows that, the integration con- 
stant M is the gravitational mass of the object provided spacetime is AF. It is 
clear that the basic requirement for an AFST is not fulfilled if A ^ because 
then one would have [5] 

e v = e- x = 1 - 2M/R - AR 2 /'S (31) 

Accordingly, in this context, Tolman dropped the A term (see p. 230) [5] and 
used (quasi) Cartesian coordinates with x 1 = x, x 2 = y and x 3 = z. Then he 
showed that the if the volume integration over B would be converted into a 
surface integral at spatial infinity, one would obtain 



/ 



B dxdydz = ip (32) 



so that 



P = J A dxdydz + ^P + J C dxdydz (33) 
By transposing and then multiplying by 2, one would obtain 

P Q = 2 J A dxdydz + 2 J C dxdydz (34) 
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Further, if the system is static, one has C = 0. In such a case one would obtain 
Po = Jv^g (r ° - Tl - If - T 3 ) d 3 x (35) 

even if the system is not spherically symmetric. Interestingly, for an AFST, 
Landau & Lifshitz[7] (pp.348, Eq.[100.19]) obtained exactly the same relation 
for the total matter plus field of a static system without invoking any pseudo 
tensor at all. This shows the physical correctness of the Einstein EMC provided 
one would evaluate it in quasi-Cartesian coordinates. 

One can see from such exercises by Tolman[4,5] and Landau & Lifshitz[7] 
that for static case, in the absence of A, one may distinctly split P into a 
matter contribution 

P (Matter) = J J=g~ T ° d 3 x, (36) 

and a pure field contribution 

Po(field) = J^gt° d 3 x = J J—g {-Tl - T 2 2 - T 3 ) d 3 x (37) 

where the integrations now effectively span only the region occupied by the 
matter. This happens because, in the absence of A, by definition, T\ = 
exterior to the body. Further, if the matter is represented by a perfect fluid, 
we will have 

T£ = (p + p)u a u b -p8 a b (38) 

where u a is the matter 4-velocity and p is the isotropic (matter) pressure. In 
the comoving frame, one obtains 

u a = 5% (39) 
and the non- vanishing components of T b a are (G = c = 1) 

T ° = p; Tl = Tl = T 3 = -p (40) 
Accordingly, in the comoving frame one has 

P (Matter) = J J~^g~ p d 3 x (41) 

and 

P Q ( field) = J 4 d 3 x = J 3p d 3 x (42) 

Much later, Herrera et al.[17, 18] used Tolman's formula for finding the 
active gravitational mass of a spherically symmetrical collapsing fluid in AFST. 
This was probably the maiden application of Tolman's generic prescription for 
a non-static system having a boundary. On the other hand, we will invoke this 
formalism to find the energy of a non-static system not having any boundary 
at all. 
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2.2 Inclusion of A For Isolated Objects? 

If the A term would be shifted to the right hand side of Einstein Eq., it follows 
that, in the comoving frame, one may replace p and p by their respective 
effective values: 

Pe = P + A/8-k; p e =p- A/8ir (43) 

Then one might think that, in such a case, the definition of the Einstein mass 
for a static isolated system would be modified into 



Po 



= J V^9(Pc + 3Pe) d 3 x (44) 



But it may not be so because of the following reasons. As seen by the exte- 
rior metric, when A is present, in a strict sense, the spacetime ceases to be 
asymptotically flat. Since, now, there would be no strict Kepler's 3rd law, one 
cannot identify the integration constant M with Schwarzschild mass. However 
it might still be possible to define "mass" by using more advanced and specific 
ideas which is beyond the scope of this discussion [13]. 

Nevertheless, it may be noted that the spacetime appears to acquire, 
strange and rather unphysical properties. Because it would be seen that while 
goo initially monotonically increases with R, it would start decreasing subse- 
quently And at R = R* , defined by the roots of the following Eq. 

200 = 1-^-^ = 0, (45) 

-ft* o 

one would have a singular situation with goo = 0! And following this, goo 
would reverse its sign to approach —00!! Since the gravitational red/blue shift 
of the photons emitted by the body depends on the value of the goo(R) at 
the point of observation, one would see strange and abnormal red/blue shifts. 
At this juncture, one would argue that, for sufficiently, large R, the metric 
of the isolated object must merge with cosmological metric, which under, the 
assumption of homogeneous and isotropy does not admit any gravitational 
red/blue shift. But the fact remains that, any mass determination might get 
swamped by a cosmological vacuum contribution 

Po(vacuum) ~ J ^g (-A/An) d 3 x (46) 

which would diverge if the spatial section would be of infinite extent. Irre- 
spective of the precise boundary condition, the foregoing Eq. indicates that, 
a completely vacuum universe should posses a net negative energy density be- 
cause of the negative contribution due to p e = —A/8n\ We would seek a more 
precise answer to this question later. 
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3 Negative Self-gravitational Energy 

If the metric determinant of the spatial section is — h, one can write 

9 = 9oo{-h) (47) 

and accordingly, following Landau-Lifshitz derivation of M for an isolated 
static system, we may write 

Pa = M = J + 3p) dV (48) 

where 

dV = V^h d 3 x (49) 

is proper 3- volume element. In modern terminology, Po — M is the ADM mass 
measured by a faraway inertial observer S°° . Thus the comoving (local) Active 
Gravitational Mass Density (AGMD) 

p g =p + 3p (50) 

indeed appears to increase due to the "weight" of pressure, It is however impor- 
tant to note that this pressure contribution is actually due to the field energy 
contribution (when computed in quasi-Cartesian coordinates): 3p = t® and 
the field energy density is positive as long as p is positive. This is similar to 
the fact that while electrostatic interaction energy can be positive or negative 
depending on the sign of the electric charges, the field energy density E 2 /8n 
is always positive [7]. Since gravitational waves carry positive energy [20-23], in 
general, we expect t{j to be positive. 

So the question which might arise is that then how one obtains negative 
self-gravitational energy? Note that the distant inertial observer S°°, however, 
sees an effective total mass density 

Pf = V9o~o~(p + 3p) < Pg (51) 

because goo < 1- And it is the difference between p g and p^ which gives rise 
to negative self-gravitational energy when viewed globally [24]. This question 
may be probed by noting that in the zeroth order post Newtonian limit, the 
Newtonian gravitational potential is obtained from 

4> = ^(.9oo - 1) « \ (52) 

Since the potential ip is negative, its coupling to the rest of the matter gives 
rise to global negative self- gravitational energy. Further, it has been explicitly 
been shown by Tolman (pp. 248-250) [5], that 



Pa = M 



J P dV + J^pipdV = JpdV-3jpdV = J{p- 3p)dV (53) 
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Thus, in this limit 

E 9~- J ZpdV (54) 

and due to global gravitational coupling, the effective AGMD seen by 5°°: 

p™*p-3p<p (55) 

Hence one has negative self-gravitational energy only when goo < 1 and in 
general goo = 9oo(R)- In such a case, as light propagates through a medium 
having negative self-gravitational energy, one would see gravitational redshift. 
But, in case, one would fix g o = 1, i.e., synchronize all clocks in a rather 
Newtonian fashion, one would not have any global self-gravitational energy. 
In such a case, neither would one see any gravitational red-shift quite like 
the Newtonian case. All red-shifts then must be due to kinematical/Doppler 
origin. Cosmological redshifts, in the paradigm of the FRW model, are indeed 
believed to be of purely kinematical origin because there is no gravitational 
redshift for g o = 1 = fixed. 

For later requirement, let us ponder, is it possible that Po — ► ever? As 
hinted by Eq.(48), in principle, it is possible in the singular limit goo — * 0. 
Only in such a case, in view of Eqs.(53) and (55), one may conceive that 
positive energy of matter/radiation/heat gets nullified by the negative self- 
gravitational energy even though field energy density t§ remains positive. Note, 
the field energy t° defined locally via Eq.(6) must not be confused with self- 
gravitational energy E g . For the latter there is no proper local definition. Thus 
even if such a vanishing of P would take place, p and p would still retain their 
original sign and there cannot be any cancellation of matter energy with field 
energy. In contrast, (negative) self-gravitational energy is essentially a global 
concept and cannot be defined at a spacetime point. Further, it is mediated by 
.9oo = goo(R) < 1- And when g o = 1, one should have E g — and then there 
would be no question of cancellation of matter energy. Thus, in the ansence of 
a variable goo < 1, an occurrence of Po = should signify p = 0. 

4 FRW Metric in Cartesian Coordinates 

The isotropic form of the FRW metric is[5](pp. 338) [8]: 

ds 2 = dt 2 - + S 3/4)2 [ rfr2 + r2 ( d02 + sin ' 0d ^ ( 56 ) 

where k can assume values of 0, +1, or — 1. This form can easily be written in 
terms of Cartesian coordinates: 

ds 2 = dt 2 - {1 f® /4)2 (dx 2 + dy 2 + dz 2 ) (57) 

where 

r 2 = x 2 +y 2 + z 2 (58) 
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Thus, 

S 2 

911 = 922 = .933 = ~J2> 5oo = 1 (59) 

and 

g ll = g 22 = g 33 = _f. g 00 = 1 {m) 

where 

f(r) = 1 + kr 2 /4 (61) 

At the beginning, let us mention that greek letters would represent x, y, z and 
in this section, repetetion of a greek letter in any expression will not imply any 
summation. Then it may be noted that 

d a f = \kx a (62) 

while 

dof = (63) 
Because of the symmetry we also have 

d a 9f3fi =p a = (64) 



4.1 Useful Christoffel Symbols 

In this quasi Cartesian basis, one may verify that since goo = 1, the non- 
vanishing connection coefficients are 



where 



Next. 



Further, 



Also, 



r a a = lg n q = ^ (65) 



q = d gu = 2 f' S (66) 



r° a , a = 4 q= Y (67) 



1 hr a 

r: a = -~ (68) 



2 / 



rl 2 = n s = -\ g " P i = \j 

All other connection coefficients vanish. 



1 Note, in the published version, there were typos in Eqs.(65) and (67) 
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5 Direct Computation of P for FRW Universe 

Let us first compute the B-term: 

5.1 Evaluation of B 

From Eq.(9), it is seen that 

®L _ r i , 1 r d , 1 r d _ r i , r ld _ p2 . r 3 _ ^ 
^TT — ~ L n + 2 ld 2 ~~ _ ~~ 12 13 _ ~/~ 

On the other hand, we have 



(70) 



dL 1 _ dL 1 _ dL 

Then one finds that 
ab dL n { dL dL dL \ dL n kx —kxf 

In general, we can write 

^ ab dL _ -kx a f 



By differentiating the foregoing Eq., we obtain, 



(73) 



By adding the three components of the foregoing equations and also by 
recalling Eq.(26), we obtain 

B = ^Uf+W) (75) 



16ttS 2 V 2 

Now plugging in the value of / from Eq.(61) in this equation, we finally find 

Further since the spatial curvature of the spacetime is 

K = ± (77) 



we may also write B as 

-K / 5 
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5.2 Evaluation of C 

From Eq.(9), one can see that when a = b, one has 



^ = ~Ka + *m ( no sum ° ver a ) ( 79 ) 



from which one obtains 



dL =n 1+ n 2 +r* m =y\=4 (so) 



and 



Similarly 



Further note that 



and 



^5ff-*01-r*02-r*03-2» « - o- 



' VI rO =+ I 5 = Z^ (81) 



c^ 1 11 2 y ^ 



2 



.22* , A _ as* (dL\_S & 



16tt 



do hriKT + S^ 1 ^, 



a.g V ^ W 



(82) 



*(w)- 3 (i-£ (83) 



5 H^J =_ ^ = S + ^ (84) 



Similarly, 

Since g 11 = g 22 = g 33 , from Eq.(27), we can write 



(86) 



By using Eqs. (83-85) in the foregoing Eq., we find that eventually C gets 
simplified as 

c 'll < 87 » 

Now let us recall the evolution Eq. of the FRW universe 

s = ~T p * (88) 

where 

p* = ( Pe + 3 Pe ) = (p + 3p)-^ (89) 
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Eqn. (88) shows that, even for the FRW case for which neither AFST nor any 
other boundary condition has been assumed, p* appears to be the AGMD. 
Then by using Eqs.(23) (87) and (88), we can reexpress 

C = ^ = § (90) 
Clearly, C has a distinct physical and geometrical meaning in terms of and 

Rq. 



6 Final form of Einstein Energy Complex 

Since for the metric (57), 



o3 

^=73 (91) 



and the FRW metric is actually a comoving metric, from Eqs. (23) and (25), 
we obtain 

By using Eq.(88), we may also write 

Clarly, one would have A = for a static case. It is seen that like the static 
case, for this dynamic case too A has a clean physical and geometrical meaning 
in terms of and Rq. The additional factor S 3 /f 3 is the ratio of the proper 
volume element to the coordinate volume element and which becomes unity for 
k = in the static case. In view of this assumed dynamic metric, eventually, 
A behaves quite unlike the case of an isolated object in an aymptotically flat 
space time where in the static case / A d 3 x = P /2[4,5]. After combining all 
the contributions, the final form of Einstein energy complex Pq for the FRW 
metric, when evaluated directly, becomes 

,0 -k / 5, 2 \ 3S/ ( S 3 \ -k / 5, 2 \ p* p« S 3 

(94) 

In particular, note that, when k = and / = 1, one has 



' : '-Al-S 3 ) = —p*(l-S 3 ) (95) 



Due to spherical symmetry, we may take the coordinate volume element to be 
dV = Airr 2 dr so that the form of Einstein energy of the FRW metric becomes 

P Q = J Anr 2 9° dr (96) 
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By inspecting, the metric (56), it becomes clear that the range of r = 0, oo for 
k = 0, 1 while r — 0, 2 for k = — 1. Further, the proper volume element is 

= rfF (97) 

Now using Eqs. (94-97) we may simplify P as 

p »=i^/( 3+ H di,+ T< v -"» w 

Here the term p»V looks like some sort of effective proper energy content. The 
1st term of the right hand side of the foregoing Eq. shows the contribution due 
to curvature of the spatial section alone while the 2nd term denotes contribu- 
tion due to spacetime curvature in general. In particular, the 2nd term implies 
a contribution with respect to a flat spacetime having no energy. Thus Einstein 
energy has a nice physical and geometrical interpretation even in this dynamic 
case. Such a physical and geometrical interpretation of the energy of the FRW 
metric has not been revealed before. In a trivial case of a flat spacetime with 
k = and V = V, obviously P = 0. 

7 Analysis of this general formula 

7.1 Position Dependence of Net Energy Density 

The first apparently anomalous thing about the Einstein energy complex (94) 
is that, it is position dependent unless k = 0. One may try to explain away this 
position dependence by telling that "gravitational field energy is not localizable 
due to principle of equivalence" . But all that this latter statement means is 
that at a given spacetime point, the value of t{j would vary as one would employ 
different coordinate systems and in Riemann normal coordinate system, one 
would have t$ = 0. It is in this sense that, one may not demand conservation 
of energy momentum in a covariant way. But once we have decided to work in 
a given coordinate system, which in this case is the quasi-Caresian coordinate 
system, and not to compare the values of t(j obtained in various coordinate 
systems, we have already compromised for such non-covariance. Such non- 
covariance however does not at all imply that a quantity which is expected to 
be spatially uniform can turn non-uniform. It is true that one speaks of non- 
localization of position of an electron/photon in quantum mechanics because 
of inherent wave-particle duality. But we are not treading into any quantum 
mechanics here. Hence we do not expect 9$ to have any inherent spatial spread 
. Consequently, in a supposed isotropic and homogeneous universe, one should 
not expect any position dependence of 9% irrespective of its the precise physical 
significance. 

Further, the Einstein pseudo-tensor need not be be blamed for such an 
occurrence which seems to defy the assumption of homogeneity and isotropy. 
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The real reason for such unexpected position dependence is the position depen- 
dence of the metric coefficients g aa . This position dependence vanishes only 
for k — when 6$ too becomes position independent. The situation here could 
be something like the following: If one would consider a spherically symmetric 
finite fluid of constant p, the pressure still would show position dependence. 
This position dependence of p can be eliminated by either (i) pushing the 
boundary to oo, or (ii) by setting p = by hand, or (iii) by pursuing the 
spatially flat Newtonian limit p — > 0. This would be clear as we would discuss 
below the case of a supposed static universe. 



7.2 Static Universe 



In the static case, the expression for Einstein energy complex becomes 



k (~ 5, 2 



0°n = 



It is known that, by means of the following coordinate transformation 

r 



1 + kr 2 /4 



the metric (56) would transform as 

df 2 

ds 2 = dt z 



+ f 2 (d6 2 + sin 2 9dtf ) 



1 — kr 2 

For simplification, we now drop bar from f to rewrite the above Eq. as 

dr 2 



ds 2 = dt z 



1 — kr 2 



+ r^(dT +sm z 6dtf) 



(99) 



(100) 



(101) 



(102) 



If we recall that the static interior Schwarzschild solution for a constant density 
may be expressed as[5,25] 



ds 2 = e v dt 2 - S 2 



dr 2 



1 - kr 2 



r 2 {d9 2 + sin 2 6dcf) 



(103) 



where k = +1 corresponds to an effective p e > 0, k = corresponds to p e = 
and k = — 1 corresponds to p e < 0, it might appear that, the FRW metric is 
similar to a special interior Schwarzschild solution where, somehow, one has 

e v = 1 (104) 

In fact, for exploring a static FRW universe, Tolman indeed arrived at Eq.(102) 
by starting from Schwarzschild interior solution (103) (see pp. 333-337) [5]. He 
found that the condition for having e v = 1 is either 



p e + 3p e = 



(105) 
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or 

Pe+Pe=0 (106) 

Essentially, Tolman probed the condition for e v to become position indepen- 
dent. On the other hand, we have probed the question from a broader perspec- 
tive where all physically meaningful quantities such as the the invariant/scalar 
acceleration experienced by an interior fluid element [25] 



a = \/-a,ia l (107) 

where 

a 1 = u^V-jU 1 (108) 

is the 4-acceleration, and radial pressure gradient p' should indeed be posi- 
tion independent. We found the following relations which, in general, suggest 
position dependence of a and p'\ 

a kSr 

(109) 



p e + 3p e y/l - kr 2 

and 



(HO) 



(Pe +Pe)(Pe + 'iPe) I - kr 2 

In a homogeneous and isotropic universe, the left hand side of the two foregoing 
equations must be independent of r. But the inspection of the right hand side 
of the same equations tells that, it could be so only when k = 0. For a static 
universe, occurrence of k = corresponds to p e = 0. It was independently 
shown that in case one would assume k = 1 , the condition that a and p' would 
not blow up at r = 1 demands that both Eqs.(105) and (106) are satisfied[25]. 
This again leads to p e = 0. Under weak energy condition, this would mean 
that a static unverse should have A = as well as mean matter density p = 0. 
Incidentally, occurrence of a mean p = need not imply absence of matter 
in cosmological models. For example, one has mean zero matter density in an 
infinitely hierarchial fractal universe[26]. But for isotropy and homogeneous, 
case, a mean p = may occur only when matter patches of the universe are 
separated by infinite distances from one another[27,28]. In any case it appears 
that in order that 9® is position independent as per the assumption of isotropy 
and homogeneity, only the k = FRW model is physically admissible. Note 
this argument is irrespective of the question whether P is conserved or not. 

8 Empty Universe with a Positive A 

As is known, for an empty universe FRW solution becomes de-Sitter solution 
ds 2 = dt 2 - S 2 {dx 2 + dy 2 + dz 2 ) (111) 

having k = and 

S = e m (112) 
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where 



H = yJA/i 

Since p = p = 0, in this case, Einstein energy complex becomes 



(113) 




) 



(114) 



Using the fact that now 



S 
S 



A 



(115) 



3 



we obtain 




) 



(116) 



But expeept for the instant of t = 0, one would have S > 1 and accordingly, 
0g < 0. Then, one finds 



Again occurrence of a negative 0q and P may be related to the fact Active 
Gravitational Mass Density of pure vacuum p e + 3j>3 is always negative for a 
positive A. However, from the concept of positivity of mass-energy, occurrence 
of a negative Po (and that too — oo) is unphysical, and this would again demand 
that A = 0! 

Further, even if one would accept the apparent result that Pq = — oo, it 
may be noted that there is latent temporal dependence of Po if A ^ 0: 



Such an occurrence would again suggest that A = if wc would expect Po to 
be conserved. 

There could be yet another reason for which one would expect A = in 
the context of de-Sitter solutions. With suitable coordinate transformations, 
the metric, can be transformed into a static form (see pp. 346) [5]: 



ds 2 = (1 - R 2 /S 2 )dT 2 - (1 - tf/S^dR 2 - R 2 (d9 2 + sin 2 Odtf) (119) 



In fact, de-Sitter originally obtained his solution in the foregoing form. As such, 
there is nothing anomalous in the fact that a non-static metric may look static 
after a coordinate transformation which involves time. But such form changes 
should not induce physical changes. Note in the foregoing metric, the radial 
coordinate R is the circumference coordinate which is a scalar and directly 
related to luminosity distances. Since now 500 = 9oo(R), one would expect 
that moving photons would experience gravitational redshift. But as per the 
FRW metric, there should not be any gravitational red-shift! Such a physical 
contradiction can be removed by realizing that A — and de-Sitter metric is 
actually flat Minkowski metric. 




(117) 




(118) 
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9 Time Dependence of Einstein Energy 

We have already discussed the issue of time-dependence of 6*{j and Po for fc = 
model in the context of de Sitter case. For more specific discussion, let us first 
note that for fc = 1, / = 1 + r 2 /4 we will have 

/■oo c3 /■oo J2 

V = 4t/ o p r2dr = ^ S *l (l + r=»/4)3 dr = 27r2s3 ' (120) 
where V is proper 3-volume. On the other hand, the coordinate volume is 

^ = oo (121) 

Also, 

P OO 

7 1 =4 7 r/ r 2 [3 + (5/4)r 2 ]dr = oo (122) 
Jo 

On the other hand, for fc = —1, / = 1 — r 2 /4 and 

V = 4tt / — r 2 dr = 4^S 3 / Tjj^dr = S 3 oo = oo, (123) 

Jo P Jo (1 - r 2 /4) 3 



and 



Then we find that 



V = [ 4irr 2 dr = ^ (124) 
Jo 

h= ( 47rr 2 [3 - {b/A)r 2 ]dr = (125) 
Jo 



Po = 0+^|(l-S 3 )(oo) = -(p*/2)(l-S 3 )(oo); k = 0, (126) 



Po = 0+-^-f (32tt/3-S 3 )(oo) = (-p,/2)(32V3-S 3 )oo; fc = -1 (127) 
and 



Po = ~ + ~{oo - 2tt 2 5 3 ) = - ~ - (p,/2)(oo - 27r 2 ^ 3 ); fc = +1 

(128) 

In case, one would consider, S = finite for the fc = 1 case, one would obtain 

Po = — oo — ^^(oo); k = +1; S = finite (129) 

In this case, Po could be +oo or — oo depending on the sign and value of p*. 
On the other hand, if S = oo for the above case, one may have 

Po = (130) 
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Since, from symmetry conditions, it has been argued that one should have, 
Po = for a closed universe [7, 29], we find that, the suposed closed case should 
actually be an open case with S — oo. In general, Po would be time dependent if 
p* would be so. For k = —1 case too, Po could be either oo or — oo unless p* = 0. 
Same is true for the simplest k = case. We have found that for an isolated 
static object in an AFST, P yields a result (in quasi-Cartesian coordinates) 
which is actually obtainable without invoking psedotensor at all. Thus, if we 
would consider that Po is a physically meaningful quantity irrespective of the 
question whether it is conserved or not and it cannot fluctuate wildly between 
+oo and — oo (even for a supposed closed geometry), we should adopt 

S/S = 0; p* = 0; S = oo (131) 

There may be a simple mathematical reason as to why one should have S = 
fixed. It is known that, if indeed S = S(t), the spacetime may be extrapolated 
back to a geometrical point. But what could be the proper 3- volume of a point: 
V(point) =?. Since this is expected to be V = 0, a time dependent open FRW 
model is immediately ruled out because in the latter case, V = oo! And of 
course, if we would demand that Po is conseved, the foregoing constaints would 
be naturally imposed on the FRW model. 



10 Discussions 

For the first time, we obtained an expression for the Einstein energy complex 
6>o for the FRW metric in a direct fashion i.e., without resorting to any super- 
potential. This direct treatment which explicitly involved T% yielded a direct 
interpretation of Po in terms of spacetime curvature and AGMD. The non- 
covaiant nature of 9® was handled by working it out in the quasi-Cartesian 
coordinates and by making no demand that its value would be same in other 
coordinate systems. Irrespective of the precise physical interpretation of 6®, 
for an assumed homogeneous and isotropic spacetime, one would expect it to 
be coordinate independent because no wave-particle like duality is involved 
here. We pointed out that such position dependence actually appears from 
the coordinate dependence of g aa (for k ^ 0). For a model static universe, 
physical obersables like "acceleration due to gravity" a and pressure gradient 
too would have similar latent coordinate dependence for k ^ 0. Thus, we find 
that though mathematically, constant curvature homogeneous and isotropic 
spacetime could be of 3 variety, physically, only the flat k = case is admis- 
sible. 

By considering the empty de-Sitter model, we found that 0§, while coor- 
dinate independent, is negative. Further P = — oo unless A = 0. Also, we 
noted that, while in one version of de-Sitter metric, one may see gravitational 
red-shift, in another version, there would be no such gravitational redshift. 
From all such considerations, we got strong hint that Einstein was right in 
rejecting A. Such a rejection need not be inconsistent with the observations 
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because there are some analysis which claims that the observations of distant 
Type la supernovae are actually in agreement with A = picture[30,31]. 

We found that for a supposed static universe, in order that, P is non- 
negative and non- infinite for the closed k = +1 case, one should have S = oo. 
Such a value of S actually tantamount to k — case because in both the cases 
curvature K = k/S 2 = 0. This again suggests that A = 0. Let us highlight a 
simple point which might have been overlooked previously. In the Rimcnnian 
normal coordinates, where d c g a b = 0, one expects all components of t b a to 
vanish. However as seen from Eq.(6) and as earlier noted by Tolman, in this 
local Lorentz frame (see Eq.[87.15])[5] 

t'i = (132) 

would not vanish unless A = 0. Thus a finite A seems to be inconsistent with 
the principle of equivalence. 

We found that 8® contains a dynamic S/S term; 

e° = e° (s,s) (133) 

In one sense, this was a welcome feature because otherwise a likely 

e° = e° (S) (134) 

alone would yield the same 8® for both a static as well as dynamic metric 
and it might appear that motion has no contribution to the total energy. 
However, it turns out that occurrence of this S/S term creates problem for 
the conservation of Pq- In such a case, P$ may vary from +oo to — oo. It 
follows then that, in order for Pq to be a physically significant quantity, one 
should restrict p* = and S = oo. Incidentally, there are claims that both 
the distant Type la supernovae and Gamma Ray Burst observations may be 
consistent with a static universe [32, 33] . In fact, it might be possible to conceive 
of static yet continually evolving universe where the patches of matter are 
infinitely separated from one another [27, 28]. Note, when p* = A = and 
p > 0, we have mean p = 0. However, such a result may be seen to be at 
variance with the patch of the universe observed now and attendant popular 
interpretations. But let us pose the question whether the observed universe is 
really isotropic and homogeneous at scales say ~ 100 Mpc? Well, Sloan Digital 
Sky Survey shows that there are structures of extent <~ 500 Mpc (Sloan Great 
Wall). Further Wilkinson Microwave Anisotropy Probe (WMAP) has found 
"void" of the extent of 280 Mpc[27,28,34]. Thus the observed universe need 
not be really described by the idealized FRW metric, and one need not be 
unduly perturbed over the suggestions gleaned from this study. In fact, cosmic 
background radition may be explained as the redshifted thermal radiation 
from eternally collapsing objects, i.e., the supposed black hole candidates [27, 
28]. However, if one would insist that the observed patch of the universe, in any 
case, must be described by FRW metric, then one might adopt two views: (i) 
Einstein's cannonical pseudo tensor does not provide a physically meaningfully 
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answer to energy momentum conservation. Such a statement would however 
overlook the fact that when evaluated in quasi-Cartesian coordinates, Einstein 
EMC does indeed find energy momentum flux for cylindrical gravitational 
waves [11]. Further, for an isolated static object in AFST, the expression for Pq 
derived from Einstein EMC indeed matches with the same derived by Landau 
Lifshitz without using any pseudo tensor at all. 

(ii) One might also adopt the view that energy conservation is invalid in 
the cosmological context. Such a view might degenerate into free violational 
of energy conservation principle in astrophysics. For instance, in such a case, 
while one might see huge eruption of energy or in some cases, one may as well 
see sudden disappearance of energy. Eruption of energy is indeed seen in as- 
trophysics in the form of supernovae and gamma ray burst and astrophysicists 
try to understand them using physics which, in the background, honour, en- 
ergy momentum consevation. And as far as sudden loss or vanishing of energy 
is considered, no such event has been recorded. In the absence of an energy 
conservation principle, motion of cosmic matter too could be unexpected and 
unpredictable. 

Thus, if we would not abandon this principle of energy conservation, there 
may be examples, how mathematically allowed dynamic motion can be forbid- 
den in GR. For instance, mathematically, one can always formulate equations 
describing contraction/expansion of clouds without any heat /radiation trans- 
fer. However, an energy conservation principle which involves global grav- 
itational energy dictates that there is no contraction or expansion of self- 
gravitating objects without heat/radiation transport [24, 35]. This principle 
however may be avoided by assuming that pressure of the fluid p = 0, in 
which case thermodynamics ceases to work. But a strict p = equation of 
state is possible only where p = - -thus such strict p = solutions would 
correspond to a fluid mass energy Pq = and the mathematical collapse would 
be devoid of any physical reality. The symmetry of the FRW metric precludes 
any heat /radiation flow and this might be forbidding motions even though 
mathematically the metric suggests motion. 

In fact, much earlier to such studies, while studying the simplified prob- 
lem of the adiabatic collapse/expansion of an uniform density sphere, Taub 
found that there would not be any collapse/expansion if the fluid would have 
an EOS [36]. Later Mansouri studied this problem with greater clarity in a 
paper entitled On the non-existence of time dependent fluid spheres in general 
relativity obeying an equation of state[37]. And the basic reason which forbade 
such motion was principle of energy conservation. Thus many general possi- 
bilities which may be initially suggested by mathematical equations, may not 
eventually be allowed by rigourous physics. Let us recall that, in general, GR 
does not allow global synchronization of clocks in presence of gravity, and on 
the other hand, clocks are expected to slow down in the presence of gravity. 
Even in the absence of gravity, special relativity does not allow global synchro- 
nization of clocks if motion is present. But in the FRW model, one expects 
to have a Newtonian like global time despite the presence of gravity and mo- 
tion! Does Einstein energy tell that such a global time is actually possible 
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only when the model would degenerate into a Newtonian one in the p — > 
limit? An honest introspection may show that, probably there were already 
strong hints that FRW metric, despite its apparent general nature, might be 
inherently, Newtonian with a mean p — > 0. This is so because, long back, it 
was shown by Milne [38] and McCrea and Milne [39] that the key results of the 
Fricdmann model can be exactly obtained by using pure Newtonian gravity. 
Nobody has ever explained why it is so! In particular, let us consider the New- 
tonian equation of motion of a test particle (not acted upon by any pressure 
gradient or other forces) of mass m lying on the surface of a uniform density 
sphere of mass M: 

F=-^ (135) 



Since 



one obtains 



S 2 

4-7T 

M = -jpS 3 , (136) 
F = mS (137) 



Since in Newtonian gravity, p ~ p/c 2 <C p, the 3p term does not appear in the 
definition of M. But if one is keen, one can introduce p here by an appeal to 
special relativity [40]. One might also push A within p. In any case, one can 
see that the FRW evolution Eq.(88) is nothing but the Newtonian evolution 
equation! In particular, the 4tt/3 factor in Eqs.(88), (135) and (138) may be 
signalling an inherent k — geometry. But, we know that GR reduces to 
exact Newtonian limit only when p — > 0. Further when gravity vanishes, the 
model must be static as found here by demanding conservation of Pq. Even 
if the cosmolgical redshift in the local patch of the universe would indeed be 
of dynamical nature, it would be a local phenomenon in such a case because 
conservation of Pq would mean an infinite static universe where such expanding 
patches are separated by infinite voids. So let us ponder whether Einstein 
energy of the FRW model is eventually explaining the mystery — why the key 
equations of the FRW model are obtainable from purely Newtonian physics? 
It is again emphasized that, in any case, observed patch of the universe filled 
with "walls" "filaments" and "voids" may not be described by the most ideal 
FRW metric. 
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